figS2_data = readtable('figS2.csv');

figure('position',[55,180,450,350]);

plot(figS2_data{:,1},figS2_data{:,2},...
    'o','MarkerSize',4.5,'linew',1,...
    'color',[0.000 0.447 0.741],'MarkerFaceColor',[0.0745,0.6235,1]);
set(gca,'fontsize',13);
xlim([0,35]);
ylim([0,0.625]);
grid on;
hold all;
plot(figS2_data{:,3},figS2_data{:,4},'-','linew',2,'color',[126,47,142]/255);
plot(figS2_data{:,5},figS2_data{:,6},':k','linew',3);
legend({'experiment','simplified model','numerical model'},'location','southeast');
xlabel('Optical pumping rate {\itR}_{pump}/{\itq} [Hz]');
ylabel('$|\mathcal{C}(\omega=-\omega_\mathrm{a})|$','interpreter','latex');
